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SUMMARY 


The  recently  derived  particular  integrals  of  the  Prandtl-Busemann 
iteration  equations  make  possible  the  extension  of  the  familiar  source - 
sink  concept  to  the  solution  of  the  higher-order  iteration  equations  for 
the  subsonic  potential  flow  over  thin  sharp-nose  symmetric  two- 
dimensional  profiles.  An  explicit  expression  is  derived  for  the  second- 
order  velocity  potential  and  velocity  components  and  a  method  for 
obtaining  the  higher-order  terms  is  indicated.  The  velocity  at  the  sur¬ 
face  of  the  Kaplan  bump  is  evaluated  to  illustrate  the  method. 


INTRODUCTION 


Iteration  methods  have  been  used  extensively  during  recent  years 
for  the  calculation  of  compressible  flows.  The  Prandtl-Busemann  or 
Ackeret  iteration  process  is  the  one  most  often  applied  in  aeronautics 
and  is  the  one  considered  herein.  In  the  past  this  iteration  technique 
has  been  used  primarily  to  calculate  the  flow  past  specific  profiles, 
each  profile  presenting  a  distinct  problem  involving  a  great  amount  of 
labor  in  its  solution. 

Van  Dyke  (reference  1)  has  given  the  particular  integral  for  the 
second-order  Prandtl-Busemann  iteration  equation  and  recently  Kaplan 
(reference  2)  has  obtained  the  particular  integrals  for  both  the  second- 
order  and  third-order  equations  and  has  described  a  method  for  obtaining 
higher-order  ones  as  well.  With  the  use  of  these  particular  integrals, 
the  higher -order  solutions  for  the  flow  over  sharp-nose  symmetric  bodies 
in  two-dimensional  flow  are  easily  obtained  by  the  method  to  be  described 
herein,  and  the  calculation  of  the  surface  velocity  or  pressure  distri¬ 
bution  for  a  specific,  profile  requires  only  an  integration.  The  procedure 
to  be  described  involves  considerably  less  labor  than  formerly  required 
and  also  has  the  advantage  of  simplicity  of  treatment  since  the  familiar 
concept  of  a  source-sink  distribution  in  the  plane  of  symmetry  is 
employed  in  the  analysis. 
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A  general  expression  for  the  second-order  velocity  potential  for 
symmetric  profiles  is  developed  and,  as  an  example,  the  flow  over  the 
Kaplan  hump  is  treated  in  detail.  The  procedure  may  he  extended  to 
obtain  the  third-  and  higher-order  perturbation  velocities. 


ANALYSIS 


The  Prandtl-Busemann  iteration  equations  for  obtaining  approximate 
solutions  to  the  exact  nonlinear  equation  governing  the  two-dimensional 
flow  of  a  compressible  fluid  are  briefly  developed  as  follows 
(reference  2): 

The  differential  equation  for  the  velocity  potential  4  is 

(1  "  ^)$xx  "  2  ^2  $xy  +  ^  ~  ^)$yy  =  0  (1) 


where 

x,y  rectangular  Cartesian  coordinates  in  flow  plane 

u, v  velocity  components  along  x-  and  y-axes,  respectively 

a  local  speed  of  sound 

and  the  subscripts  x  and  y  denote  partial  differentiation  with 
respect  to  these  designated  variables.  With  the  introduction  of  a 
characteristic  length  c/2,  where  c  is  the  chord,  and  the  undisturbed 
stream  velocity  U  as  the  unit  of  velocity,  the  quantities  x,  y,  u, 

v,  and  <J>  are  used  hereinafter  to  denote,  respectively,  the  nondimen- 

sional  quantities  2x/c,  2y/c,  u/U,  v/u,  and  2$/Uc  while  a  retains 

its  original  meaning. 

Let 

Y  =  tY±(x)  +  t2Y2(x)  +  0(t3)  (-1<  x  <  l) 

and 

Y  =  0  (x  >  1)  x  g  -1) 


where  0(t3)  denotes  terms  of  the  order  t^,  define  a  slender  sym¬ 
metric  body  of  thickness  coefficient  t  lying  on  the  x-axis  between  -1 
and  1. 
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The  appropriate  boundary  conditions  for  the  flow  over  a  two- 
dimensional  body  in  an  unbounded  stream  are: 

At  inf inity 


and  on  the  body 

®y(x,Y)  =  ®x(x,Y)§  (2) 

It  is  now  assumed  that  $  may  be  expanded  in  the  form 

$  =  X  4-  0^  +  0£  4-  •  • 

where  0n+q  and  its  derivatives  are  small  compared  with  0n,  and  0n 
is  of  the  order  tn.  Then,  when  the  expansion  for  <J>  is  substituted 
in  the  differential  equation  (l),  the  linear  equations  for  0^  and  02 

are 

P201xx  +  01yy  =  0  (.3) 

0202xx  +  02yy  =  2Moo2  |u  +  a)P201xx01x  +  01xy01y]  ■  (*0 

M  ^ 

where  a  =  —  *  ^  P2  =  1  -  McxP,  Moo  is  the  Mach  number  at  infinity. 

:  .  13 

and  7  is  the  ratio  of  specific  heats  at  constant  pressure  and  constant 
volume. 

The  boundary  conditions  for  0]_  and  02  >  to  the  order  t  and  t^, 
respectively,  are: 

On  the  body 


02y(x,O) 
and  at  infinity 

01x  =  ^ly  =  ^2x  =  ^2y  =  0 


(7) 
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The  well-known  source -sink  solution  of  the  differential  equation 
for  valid  for  symmetric  profiles  and  satisfying  the  boundary  con¬ 

ditions  given  by  equations  (5)  and  (7),  is  (see  reference  3>  for  example) 


0l(x>y)  =  f  ^  Ip1  logg'ju  -  i)2  +  P2^  d? 


The  solution  of  the  nonhomogene ous  differential  equation  for  02j 
equation  (4),  can  be  expressed  as  the  sum  of  the  particular  inte-  . 
gral  t2(x,y)  and  a  function  qp2(x>y)  satisfying  the  homogeneous 
equation. 


p2cp 


2xx 


+  <P, 


2yy 


=  0 


The  particular  integral  equation  (4)  (reference  2)  is 


*2(x>y) 


=  M»20 


lx 


4 ♦ §  y. 


(9) 


With  the  use  of  equations  (8)  and  (9)>  *2x  and  Qxe  easily  shown 

to  vanish  at  infinity  and  the  boundary  conditions  for  cpg;  froD1  equa¬ 
tions  (6)  and  (7),  become: 


On  the  body 


q^XjO) 


dF(x) 

dx 


where 


M“2(1  + 


and  at  infinity 


<P 


2x 


It  is  to  be  understood  that  all  differentiations  are  made  prior  to 
setting  y  equal  to  zero. 
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The  differential  equation  and  boundary  conditions  for  92  8X6 
the  same  form  as  for  0]_.  Therefore,  the  solution  for  02  niay  be 
expressed  as 

0o(x>y)  =  -h  f  ^4plog e\|(x  -  i)2  +  p2^  d6  +  io(x,y) 


provided  that  J  |F(|)|d|  converges  and  F(i)  is  sectionally  con- 

tinuous.  Since  the  source -sink  solution  may  he  obtained  by  using  the 
concepts  of  Fourier  integrals,  it  follows  that  the  function  F(£)  must 
satisfy  the  conditions  which  are  required  for  a  function  to  be  repre¬ 
sented  in  a  Fourier  integral.  These  conditions  are  sufficient  #to  insure 
the  existence  of  the  solution  although  somewhat  broader  conditions 
for  F‘U)  are  known.  These  conditions  seem  to  exclude  the  representa¬ 
tion  of  the  second -order  flow  past  blunt  bodies  by  a  continuous  distri¬ 
bution  of  sources  and  sinks  on  the  axis  alone. 

2  ' 

Then  the  velocity  potential  to  the  order  t  is 
0(x,y)  =  x  +  J ^  t  loge\|(x  -  I)2  +  d|  + 

i2(x>y)  +  0(t3)  (10) 


and  the  velocity  components  on  the  body  to  the  same  order  are 


u  =  4>x(x,y) 


1  +  i=-  P 


dYl  ,  dF(|)  |  dl 


+  tY101xy(x>°)  + 


m„2(i  +  | 


0ixx(x>o)0l(x,°)  +  [0ix(x/OT]  \  +  0(t3) 


v  =  4>v(x,Y) 


=  t  Sr]!  +  01*(x’°ll  +  } 


2  dY2 


+  O(t^) 


where  P  denotes  the  Cauchy  principal  value,  of  the  integral. 
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Once  02.  is  determined  for  a  given  profile ,  the  second-order  ^ 

velocity  potential  and  the  corresponding  velocity  components  can  he 
found  from  equations  (10),  (ll),  and  (12). 

The  third-  and  higher-order  terms  for  the  velocity  potential  or 
velocity  components  can  he  determined  in  the  same  manner  as  the  second - 
order  terms.  The  third-order  term  for  the  velocity  potential  is  written 
as  the  sum  of  a  function  satisfying  the  homogeneous  differential  equa¬ 
tion  and  a  particular  integral  (from  reference  2).  Then  the  boundary 
conditions  for  the  homogeneous  equation  are  again  in  a  form  similar  to 
those  for  0q  and  the  third-order  solution  is  the  sum  of  the  particular 
integral  and  an  expression  for  a  source  distribution  with  the  source 
strength  determined  by  the  boundary  conditions. 

It  appears  that  this  method,  as  developed  herein,  is  not  applicable 
for  evaluating  the  flow  over  all  bodies.  For  blunt  bodies,  the  method 
must  be  refined  to  include  other  solutions  in  conjunction  with  the 
source-sink  solution  in  order  to  satisfy  the  boundary  conditions.  It 
might  be  noted  that  the  troublesome  terms  in  F(l)  arise  from  the 
particular  integral.  Since  the  particular  integral  is  multiplied  by  M^, 
the  difficulty  encountered  is  a  Mach  number  effect. 

It  does  not  seem  possible  to  obtain  a  second-order  solution  for  the 
flow  past  a  sharp- nose  body  at  a  small  angle  of  attack  (to  the  second  order 
in  thickness  and  first  order  in  angle  of  attack)  by  only  a  source-sink  and 
vortex  distribution  on  the  axis.  An  integral  equation  must  be  solved  to 
obtain  the  function  in  the  lift  problem  corresponding  to  F(|)  in  the 
thickness  problem.  Therefore,  it  is  difficult  to  make  any  general  state¬ 
ments  about  the  conditions  this  function  must  satisfy.  However,  the 
particular  integral  seems  to  introduce  a  singularity  in  the  expression 
for  the  lift  coefficient  which  cannot  be  canceled  by  the  solution  to  the 
homogeneous  equation.  Here  again  the  method  must  be  refined  in  order  to 
obtain  the  second-order  lift  solutions. 

In  general,  the  integral  for  02  or  its  derivatives  are  more 
difficult  to  evaluate  than  those  for  0q  and  in  some  cases  it  may  be 
desirable  to  approximate  the  integral  for  02.  When  0q  is  known,  the 
integrand  of  the  integral  for  02  or  its  derivatives  may  be  expanded 
in  a  Taylor  series  and  the  integration  performed  term  by  term. 

The  Taylor  expansion  of  the  function  ^  ^ )  about  the  point  x  is 

dS 

00 

Mi)  = 
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where 


tSt1  -  *8  i2  ♦  e2«i^(*,o)  ♦  32t  |i  0lx(*,o) 


2/i  ^  O 


1  +  ^Ixy^*'0^*'0) 


F^(x)  =  —  F(x) 
dx11 


Then,  when  ^  of  equation  (13)  is  substituted  in  equation  (11)  and 

the  resulting  equation  is  integrated  term  by  temij  the  u  component  of 
velocity  becomes 


1  +  0-1  Y(x,0)  +  t2Yn  ^ 


1  dx 


M„2(l  +  |)  0ix2(x^0)  +  0wY(x,O)01(x,O)  - 


L-  n=l 


(1  -  x)n  -  (-1  -  x)n 


(-1  <  X  <  1 


CALCULATION  OF  THE  FLOW  PAST  THE  KAPLAN  BUMP 


As  an  example  of  the  use  of  the  second-order  equations,  the  velocity 
at  the  surface  of  the  Kaplan  bump  is  evaluated.  The  equation  of  the 
bump  is 


Y(x)  =  t(l  -  x2)3/  -  3x2t2(l  -  x2)3/2  +  0(t3) 
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From  equation  (8),  the  expression  for  0ix(x>°)  is 

n  1 


cLYi  d|  3t 


1  l\jL  -  I2 


- 1« p  J_,  if  -  -g p  J.,  « 

By  the  substitution  of  i  =  cos  0,  the  expression  for  0lx(x,O)  becomes 

ae  • "  /0”  r^ra) 

With  sin20  replaced  by  i(l  -  cos  20)  and  with  the  use  of  the 
following  integral 


cos  n0  d0  -it  sin  na 


cos  a  -  cos  6  sin  a 


the  expressions  for  01y(x,O)  and  0^(x,O)  are 


Wx<°>  ■  (I  -  \) 


-  |(f 


X  -  X 


The  u  component  of  velocity  (equation  (11))  is 


u  =  4>„(x,Y) 


=  1  +  ■=-  P 
3* 


-1 


dYl  dFU) 
dT  ■  d| 


M„2(l  +  |)[0lxx(x,O)01(x,O)  +  +  0(t 


(14) 


(  -1  <  x  <  1) 
(  -1  <  x  <  1) 


(15) 


dF(x)  _  2 

__  (3  t 


Yi0ixx(x>°)+^r  + 1 

M*,2^  +  |)01xy(x^°)^l(x^°) 


2  ^2 

chc 


with 
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'  9 


d^Y  ( x ) 

where  01v,r(x,O)  =  t  - — .  The  various  ’  terms  arising  in  I  are 

d.  CLX 


lxy 


dx 


Yl0lxx(x,O)  =  -^(1  -  x2)3/2  =  ^(x3  -  x)\/l  -  x2 


dY 

dx 


i  ^(x,0)  .  -f(l  -  x3)^7? 


dY 


^  =  (-6x  +  15x3)\/l  -  x2 


d2Y 


-  0i(x,O) 


-3t^  x  -  l^x3  +  2x^  j 


i 


dx  \/l  .  x 

The  following  integrals  will  he  required: 
„1 
•1 


(16) 


p  i  £  Vl  -  I2  dg  _  p  Pn  cos  9  sin2S  d@ 


x  -  5 


x  -  cos  6 


'  *(-5  +  *2) 


P  r"^  e  3  Vl  -  dl  p  f  *  cos^9  sin^0  dfl 

J_  1  x  “  ^  Jo  x  "  cos  9 


1  XT  V 
TB  “  *2"+  x 


%  dl 


=  P 


cos  9  d 9 


-1  \)l  ~  f(x  -  |)  J 0  x  -  cos  9 


-  -n 


1  g3dg 


-1  )|l  -  S2(*  -  I) 


=  p 


n  cos3g  d.9 
x  -  cos  9 


■  -”(!  +  *2) 


E5dg 


= p  / 
y  0 


*  cos 5  9  d9  / 3  x2  4) 
x  -  cos  9  =  "n  H  +  T  +  X  ; 


r  (17) 


1  -  l"(x  -  |) 
where  repeated  use  has  been  made  of  equation  (14). 

From  equation  (15),  with  the  use  of  equations  (l6)  and  (17),  and 


with 


Y1  — 7^-  =  -3t2(l  -  3x2  +  2x^) 
dx2  v  7 
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2  +  01(^,o)0lxx(:x,o)  =  pig  .  &2  +  5*4) 

the  x-component  of  velocity  on  the  body,  to  the  order  t2,  is 
u(x,Y)  =  1  -  l  £(2x2  -  l)  +  ^kx4  -  3x2  +  |  +  53x4  -  |  px2  +  i  p  - 


3  *2  +  M  2(3  4  _  3  2  1 

8  p  00  a\2  2  X  +  4 

in  agreement  with  reference  4. 
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Abstract 


The  recently  derived  particular  integrals  of  the 
Prandtl-Busemann  iteration  equations  make  possible  the 
extension  of  the  familiar  source- sink  concept  to  the 
solution  of  the  higher-order  iteration  equations  for 
the  subsonic  potential  flow  over  thin  symmetric  two- 
dimensional  profiles.  An  explicit  expression  is 
derived  for  the  second-order  velocity  potential  and 
velocity  components  and  a  method  for  obtaining  the 
higher-order  terms  is  indicated.  The  velocity  at  the 
surface  of  the  Kaplan  bump  is  evaluated  to  illustrate 
the  method. 
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